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The radius of gyration and the hydrodynamic radius of a number of linear, double-stranded 
DNA molecules, ranging in length from 3 to 300 kilobase pairs, extracted from modified 
strains of E. coli and dissolved in the presence of excess salt in a commonly used solvent, 
Q \ have been obtained by static and dynamic light scattering measurements. Using scaling re- 

lations from polymer solution theory, the ^-temperature for this solvent, and the crossover 
^ I swelling of the two radii have been determined. It is shown that, as in the case of neutral 

^ I synthetic polymer solutions, dilute DNA solutions exhibit universal swelling behaviour. 

^ ■ Comparison of the universal swelling with previously obtained results of Brownian dy- 

(N 

namics simulations enables an estimate of the solvent quality of a DNA solution at any 
temperature and molecular weight. The zero shear rate viscosity of semidilute unentan- 
gled DNA solutions, measured across a range of molecular weights and concentrations at 
;h \ a fixed value of solvent quality, is shown to obey power law scaling with concentration 



predicted by scaling theory, but with an eff"ective exponent Vgff that depends on the solvent 
quality, in place of the Flory exponent v. 
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I. INTRODUCTION 



The Theological behaviour of semidilute solutions of linear uncharged polymers has not been 
as extensively characterised experimentally or theoretically as the behaviour of dilute and concen- 
trated solutions. The most successful description so far of the equilibrium static and dynamic 
properties of semidil ute polymer solutions comes from s caling theories that are based on the 



concept of a "blob" 



Rubinstein and Colb 



J 



de Gennes. 



19791 : 



Doi and EdwardsL 



1986 



Grosberg and Khokhlov 



1994 



200311 ■ Indeed results of recent investigation s of the rheological behaviou r 



2005 



200811 



of semidilute solutions have been interpreted within this framework HHeo and Larsonl . 
It should be noted, however, that power law descriptions of scaling behaviour derived from blob 
theory are expected to be valid only in the limits of 6 and very good solvents. More recently, 
the development of advanced mesoscopic simulation techniques that are capable of capturing the 



many body interactions between segments of the poly mer chains has enab 



the behaviour of s emidilute solutions for the first time [[Ahlrichs et al. 



Stoltz et al. ■ 



2001 



ed the simulation of 



Huang et al. 



2010 



2006D. The use of a repulsive pair- wise excluded volume potential between polymer 
segments in these simulations, however, implies that their results are restricted to the good solvent 
regime. The rheological behaviour of semidilute solutions in the crossover regime between 6 and 
good solvents has consequently not been systematically studied so far. In this paper, we carefully 
experimentally investigate this regime and examine in particular if the zero shear rate viscosity 
exhibits universal scaling behaviour. 

Semidilute solutions of DNA have been used to carry out the studies reported here. DNA is 
a polyelectrolyte, but under high salt condit ions behaves like a charg e neutral worm-lik e chain 



with a persistence 



Marko and Siggia . 



length of about 50 nm IIBarrat and Joanny 



1995 



Robertson et al. 



20061 : 



Smith et al.i 



1996 



1996al : 



Bustamante et al.i 



Sorlie and Pecoral . 



19941 : 



19900 



The monodispersity of DNA and it's ability to be stained for visualization has made it an excellent 



from equi 


ibrium and at a range of concentrations [iFuiimoto et al. . 


1994: 


Hodnett et al. 


1976; 


Laib et al. 


2006: MaratI 


lias et al. . 2000: Navvelt et al. . 


2007: 


Nicolai and J 


vlandel. 


198( 




^ecora. 


1991: Ross and Scruggs 


1968: Sibilevaera/.. 1987: 


Smith and Chu. 


1998: 


Valle et al. . 


2005]. It 



has consequently been used extensively for a variety of experimen t al measurements. For instance 



meas u rements have been made of the diffusivit y 



19961 : 



Langowskil 



1987 



Robertson et al. 



20061 : 



llChirico et al. 



19891: 



Sells and Pecoral . 



1995 



'ishman and 



Smith et al 



atterson. 



mm. 



intrinsic viscosit y IIDotv et al. 



dilute solutions IHur et al. 



1958 



2001 



Leighton and Rubensteinl 



Liu et al. 



196911 . rela xation times in semi- 



2009 1, extensional viscosity llSunthar a/. , 



2003 



Schroeder et al. 



200511 



200311 



and conformational evolution in cross-slot flows [[Babcock et al. . 
In all these experimental observations, the precise quality of the solvents that were used has not 
been reported. However, in order to characterise the crossover behaviour of semidilute solutions, 
it is essential to determine both the ^-temperature, and the solvent quality of the solution at any 
temperature above the ^-temperature. 

In this paper, we find the ^-temperature of a DNA solution with the help of static and dynamic 
light scattering measurements. The methodology for determining the solvent quality of the solu- 
tion is also established by systematically characterising DNA molecules in an extended range of 
molecular weights from 3-300 kilobase pairs (kbp) (« 10^-10^ daltons), and at a variety of dif- 
ferent temperatures. Basically, we show that, similar to the behaviour exhibited by charge neutral 
synthetic polymer solutions, the swelling of DNA solutions in excess salt also exhibits universal 
behaviour in the context of two properties: (i) the swelling of the radius of gyration (Rg), which is 
a static equilibrium property, and (ii) the swelling of the hydrodynamic radius (Ru), which is a dy- 
namic equilibrium property. By superimposing the experimental master curves on exact Brownian 
dynamics simulations predictions of the swelling, the solvent quality of the solution at any given 
temperature and DNA molecular weight is determined. 

With the solvent quality determined, measurements of the zero shear rate viscosity have been 
carried out for DNA of different molecular weights, and at different temperatures, ensuring that all 
solutions have identical solvent qualities. We find that the dependence of zero shear rate viscosity 
on concentration is in perfect agreement with the power law scaling predicted by blob theory in the 
limits of 6 and very good solvents. More interestingly, we show that the behaviour in the crossover 
regime can also be represented by power laws in terms of an effective exponent Veff that depends 
on the solvent quality. 

The plan of the paper is as follows. In the next section, we briefly describe the protocol for 
our experiments, with details deferred to the supporting information. The procedure used here to 
determine the solvent quality of a DNA solution is discussed in section Hill with a subsection on 
the method adopted here for determining the 6 temperature of the solution (IIII Al) . and another in 
which the swelling of DNA is compared to previously reported swelling of charge neutral synthetic 
polymers and to the results of prior predictions by Brownian dynamics simulations (|IIIBI) . In 
section |Wl the concentration dependence of the zero shear rate viscosity in dilute and semidilute 



solutions is reported, and the existence of universal crossover behaviour in the semidilute regime is 
demonstrated. A comparison of the persistence length of DNA measured here with several earlier 
observations is reported in Appendix|Al and the reliability of the current measurements under poor 
solvent conditions is discussed in Appendix |Bl Our conclusions are summarised in section IVl 



II. METHODOLOGY 



A. Light scattering 



In order to experimentally characterize the properties of double- stranded DNA, a range of large 
molecular weight DNA, each with a monodispers e population, is de sirable. The present work is 
greatly facilitated by the work of Smith's group IILaib et alx 1200611 who genetically engineered 
special double-stranded DNA fragments in the range of 3-300 kbp and incorporated them inside 
commonly used Escherichia coli (E. coli) bacterial strains. These strains can be cultured to pro- 
duce sufficient replicas of its DNA, which can be cut precisely at desired locations to extract the 
special fragments. 

The E. coli stab cultures were procured from Smith's group and the DNA frag ments were 
extracted, linearized and purified ac cording to standard molecular biology protocols IILaib et al. , 



2006 



Sambrook and Russell! 



200111 . The various DNA fragments used here are listed in Table I 



of the supplementary material. In addition to the DNA samples procured from Smith's group, 
two low molecular weight DNA samples (2.9 and 8.3 kbp), procured from Santosh Noronha's 
laboratory at IIT Bombay, have been used in the light scattering measurements. The supplementary 
material also includes details of working conditions, procedures for preparation and quantification 
of linear DNA fragments, particulars of the dynamic and static light scattering measurements, and 
the approach used here for the determination of second virial coefficients. 



For each molecular weight, the purified linear DNA pellet was disso 



EDTA Buff"er), whic h is commonly used in physic s IIRobertson et al. . 



Smith et al. 



1996all and rheology llSunthar et alx l2005ll experiments. It contains 0.5 M NaCl, 



2006; 



ved in a solvent (Tris 



Smith and Chul 



1998 



which we establish short ly below is above the threshold for observing charge- screening eff"ects 
I Marko and Siggiai 1995 1. Consequently, the DNA molecules are expected to behave identically 
to neutral molecules. The detailed composition of the solvent is given in Table II of the supple- 
mentary material. 



We have used Static Light Scattering (SLS) (BI-200SM Goniometer, Brookhaven Instruments 

Corporation, USA) to determine Rg and the second virial coefficient Aj at diff'erent temperatures 

(from 10-35°C), and dynamic Light Scattering (DLS) (Zetasizer Nano ZS, Malvern, UK) has 

been employed to determine at difl^erent temperatures (from 5-35°C). A concentration of c/c* 

= 0.1 (for DLS) and c/c* =0.05-0.1 (for SLS) has been used throughout this study. Note that, 
M 

c* = 7, is the overlap concentration, calculated based on the measured values of Rg, 



NAi4n/3)Rl 



with Na being the Avogadro number. Typical c* values for different DNA molecular weights are 
listed in Tabl e III of the s upple mentary material. Initial estimates of c* were based on values of Rg 



suggested in 



Laib et al. 



[I2OO6I1 for molecular weights that are identical to those used in this study. 



B. Rheometry 

The viscosity measurements reported here have been carried out on three different DNA molec- 
ular weight samples, (i) 25 kbp, procured from Smith's group as described above, (ii) linear ge- 
nomic DNA of /1-phage (size 48.5 kbp), purchased from New England Biolabs, U.K. (#N301 IL), 
and (iii) linear genomic DNA of T4 phage (size 165.6 kbp), purchased from Nippon Gene, Japan 
(#314-03973). 

All measurements were carried out with a Contraves Low Shear 30 with cup and bob geometry 
(IT/IT). This instrument is efficient for measuri ng low viscosities and h as very low zero-shear rate 



viscosity sensitivity at a shear rate of 0.017 s llHeo and Larsonl 



of this device is detailed in an earlier study [Heo and Larsonl 



200511 . The measuring principle 



200511 . The primary advantage of 



using this is a small sample requirement (minimum 0.8 ml), which is ideal for measuring DNA 
solutions. The zero error was adjusted prior to each measurement. The instrument was calibrated 
with appropriate Newtonian Standards with known viscosities (around 10, 100 and 1000 mPa-s 
at 20°C) before measuring actual DNA samples. Values obtained fall within 5% of the company 
specified values. 

Steady state shear viscosities were measured at a temperature range of 10 to 35°C for all lin- 
ear DNA samples and a continuous shear ramp was avoided. Prior to measurements, A and T4 
DNA were kept at 65°C for 10 minutes and immediately put in ice for 10 minut es for their maxi 



murn concentrations. This was done to prevent aggregation of long DNA chains [|Heo and Larson . 



200511 . The shear rate range of the instrument under the applied geometry is from 0.01 to 100 s 



-1 



At each shear rate, a delay of 30 seconds was employed so that the DNA chains get ample time to 



relax to their equilibrium state. At each temperature, a 30 minutes incubation time was employed 
for sample equilibration. 

The shear viscosities at each temperature and concentration in the plateau (very low shear 
rate) region were least-square fitted (along with the errors) with straight lines and extrapolated to 
zero shear rate. The zero shear rate viscosities were determined in this manner for the different 
molecular weight DNA samples. 



III. DETERMINATION OF THE SOLVENT QUALITY OF A DNA SOLUTION 

In all experimental observations carried out so far on DNA solutions in excess salt, it is im- 
plicitly assumed that their behaviour is identical to the behaviour of neutral polymer solutions. 
Indeed all experimental observations so far support this implicit assumption in terms of general 
expectations of their behaviour. However, a unique aspect of solutions of synthetic polymers of 
sufficiently large molecular weight, namely, the existence of universal behaviour due to the self- 
similar character of polymer molecules, has not been investigated for DNA solutions. Universal 
behaviour enables the representation of experimental observations in terms of master plots, and 
the description of solution behaviour in terms of scaling variables. 

An example of universal behaviour in the dilute limit that is of relevance to this work is that 
exhibited by the swelling of individual polymer chains when the polymer solution temperature T 
increases above the 6 temperature, Ty, due to repulsive intramolecula r interactions. Such expansio n 



200311 ■ 



is well described by dimensionless expansion (or swelling) ratios [[Rubinstein and Colby , 
for instance, for Ro and ch for Ry^: 

Rg(T) 
g 

Ru(T) 

«H = —f- (2) 

where, Rg(T) and RniT) are the radius of gyration and the hydrodynamic radius at any temper- 
ature T > To, and i?g and i?^ represent the 0-temperature values of these quantities. As is well 
known, while the scaling of the radius of gyration and the hydrodynamic radius with molecular 
weight M obeys universal power laws in the special cases of 6, athermal, and non-solvents, in the 
intermediate crossover region between 9 and athermal solvents, the ratios Cg and Ch can be col- 
lapsed onto master plots when plotted in terms of a single scaling variable, r = ^1 - ^ j Vm, that 



combines the dependence on temperature and molecular weight HHayward and Graessleyl . 11999 



Rubinstein and Colby . 



2003h 



In the case of neutral synthetic polymer solutions, this universal behaviour is commo nly inter 



1981 



Yamakawa . 



200 ih . Accord- 



preted in terms of the two-parameter theory p Miyaki and Fujital . 
ing to the two-parameter theory, the values of properties at any temperature T under good solvent 
conditions, can be expressed as an unique function of just two parameters, namely, the radius o ' 



gyrat ion under 6 conditions, i?g, and a solvent quality parameter, commonly denoted by z HSchafei 
199911 . The parameter z is proportional to the square root of the number of monomers in the poly- 
mer chain, with the proportionality constant being related to the binary cluster integral, which 
reflects the departure of the temperature from the 0-temperature, and th e chemistrv of solvent 



mediated interact i ons b etween monomers in different parts of the chain IIRubinstein and Colby . 



2003 



Yamakawal . 



200 ih . Thus, for instance. 



R,iT)=Rlf,iz) 



(3) 



implying the universal two-parameter theory prediction, org = fg{z), where /g is a universal func- 
tion of z, whose specific form will be elaborated later. Similar predictions are made by two- 
parameter the ory for q;h and the swelling of th e viscometric radius R^, which is based on the intrin- 



sic viscosity iHayward and Graessleyl II999I1 . The connection between universal two-parameter 



theory predictions and the master plots obtained by experiments on synthetic polymer solutions 
has been made by relating z to the scaling variable t through a simple linear relation z = kr . 



where k is a co nstant that depends on polymer- solvent chemistry [|Hayward and Graessley 



Schafei 



199911 . 



19991 : 



Away from equilibrium, a scaling variable in addition to z is required to characterize dilute 
solution behaviour. For instance in shear flows, it has been found that including a characteris- 



tic shear rate ^ enab 
[ Kumar and Prakas h i 



es the representation of rheological behaviour in terms of universal plots 



2004 



Prakash . 



200211 . while in extensional fl ows, the inclusion of the Weis- 



senberg number Wi enables a similar representation [ISunthar and Prakashi . 



titative comparison of the stretch evolution of DNA in a cross slot flow llSmith and Chul 



200511. Indeed, a quan 



and the extensional viscosity of DNA solutions in uniaxial elongational flow jSunthar et ali 



998 



20051) 



was shown to be possible once an estimate of the solvent quality z of the solution was made. 
However, to the best of our knowledge, a systematic study to estimate the solvent quality of DNA 
solutions has not been attempted so far. The characterization of DNA solutions in terms of the two- 



7 




FIG. 1. Dependence of /?h on salt concentration for two different molecular weights at 25 °C. 



parameter theory would provide a powerful predictive tool for theoretical analysis of its dynamics 
in equilibrium and in flow. 

In order to proceed with the experiments proposed in this work, it is essential to first establish 
the salt concentration regime in which DNA (a polyelectrolyte) behaves as a neutral polymer. We 
have measured the hydrodynamic radius of two different linear DNA fragments across a range of 
salt concentrations (from 0.001 to 1 M) at 25°C and the results are displayed Figure [B It is clear 



from the figure that complete charge screening occurs above 10 m M NaC 



with earlier dynamic 



Soda and Wada . 



. Thi s is in agreement 



1987 



Liu et al. 



20001 : 



ight scattering studies on linear DNA [iLangowskiL 
198411 . Since the solvent used here contains 0.5 M NaCl, the light scattering 
experiments of the current study are in a regime well above the threshold for observing charge 
screening effects. We have also determined the persistence length of DNA under conditions where 
the charges are completely screened, and our estimate agrees with earlier measurements of the 
persistence length under such conditions. These results are discussed in Appendix lAl 

Figure [2] compares present measurements o f the dependenc e of hydrodynamic radius on 



molecular we ight with previous measurements [[Robertson et ali 



25°C. While 



Smith et al. 



2006 



Smith et al. . 



1996all at 



[|l996a|l used fragments and concatenates of A phage DNA to obtain 



molecules across the w ide range of molecular weights that were studied, the measurements of 



Robertson et al. 



[I2OO6I1 were carried out on molecules identical to those that have been used here. 
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FIG. 2. Comparison of the molecular weight de pendence of hydrodynamic radius, obtained previously by 



Smith fl/.i Iil996all and 



Robertson et al. 



1 200611 at 25°C with the current work. 



These results were used by Smith and co-workers to demonstrate that double- stranded DNA in so- 
lutions with excess salt obey the molecular weight scaling expected from polymer solution theory 
in the presence of excluded volume interactions. In the present work, we extend the measurements 
to a wide range of temperatures in order to carry out a systematic study of the crossover region 
between 9 and very good solvents. 

A. The ^-temperature 

The ^-temperature for a particular polymer- solvent pair can be inferred from the measured data, 
by several signatures. We have employed two methods — molecular weight scaling and the second 
virial coefficient transition. At the 0-temperature, Rg and R^ scale as for large molecular 
weight linear polymers. With this expected behaviour, two dimensionless scaling variables can be 
defined: 

n„ ^ (5) 

a yM 
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FIG. 3. Scaled hydrodynamic radius (IIh) as a function of molecular weight in base-pairs or bp. Slopes are 
calculated at each temperature by a least-squares linear fit of the last four data points. 
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FIG. 4. Scaled radius of gyration (Fig) as a function of molecular weight in bp. Slopes are calculated at each 
temperature by a least-squares linear fit of the last four data points. 
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FIG. 5. Variation of slope obtained from Figure [3] with temperature. Readings corresponding to 5 to 25°C 
are fitted with a second order polynomial to calculate the 0-temperature. 



where, a is a constant with dimensions of length, which we have set equal to 1 nm. Both IIh and 
rig should increase with molecular weight for good solvents, remain constant for theta solvents, 
and decrease for poor solvents. 

Figures |3] and |4] show IIh and Ilg plotted against M in log-log plots. In line with expectation, 
the data shows a monotonic trend: decreasing with molecular weight for certain temperatures, 
remaining relatively constant at intermediate temperatures, and increasing with M for higher tem- 
peratures. In order to determine the ^-temperature precisely, the effective slope of the curves at 
sufficiently large molecular weights was plotted against absolute temperature and fitted either with 
a second order polynomial or subjected to linear interpolation. The temperature where the sign 
of the slope changes from negative to positive, is taken to be Tq. A value of Tgi = 13 + 2 °C was 
estimated for the DNA-Solvent pair, after taking into account experimental errors in Ilg or IIh (see 
Figure [5]). 

Figures [3] and m also indicate that for the relatively smaller molecular weight samples, at all the 
observed temperatures, both IIh and Ilg are roughly constant with increasing molecular weight. 
This behaviour is consistent with the hloh scaling argument which anticipates that below a length 
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FIG. 6. Variation with temperature of the ratio of second virial coefficient to molecular weight. Data points 
are least-squares linear-fitted to calculate the 0-temperature. 



scale corresponding to a thermal blob, the energy of excluded volume interactions is less than k^T 
(where is the Boltzmann constant), as a result of which chain segments with size less than the 



Rubinstein and Colby 



blob s ize obey ideal chain statistics, and scale as M"^ at all temperatures 
1120031]. In Appendix IHl we exploit the consequences of this scaling argument to determine the 
size of thermal blobs as a function of temperature in DNA solutions, under good and poor solvent 
conditions. 

An alternative method employed in synthetic polyme r literature for est i mating the theta temper- 



ature i s by determining t he sec ond virial coefficient A2 IIKrigbaum . 



1963 



Shultz and Flory , 



1954 



Orofino and Mickey Jr. . 



195211 . either from the concentration dependence of osmotic pressure. 



light scattering from dilute polymer solutions, or by measuring intrinsic viscosities of dilute poly- 
mer solutions. A? is a direct measure of excluded volume interactions between pairs of chains 



jRubinstein and Colby . 



20031] . In this work, a quantity that is linearly proportional to the sec- 



ond virial coefficient is determined for a single, linear, medium molecula r weight DNA fragment 



(25 kb p) by static light scattering, according to the procedure described in 



Orofino and Mickey Jr. 



11 196311 . The exact values of A2 are not determined since we are interested only in the tempera- 



12 



ture where it vanish es. The sign of A2, which is i ndicative of repulsion or attraction between the 
segments of a chain IIRubinstein and Colbyi . 1200311 . is positive in good solvents, negative in a poor 
solvents and vanishes at the theta temperature. The temperature at which A2 vanishes, determined 
by interpolating the measured data, is used here as an estimate of Tg, as shown in Fig. [6l We find 
Tg = 13 + 2 °C for the solvent used in this study. This value agrees, within experimental error, 
with the value obtained using the molecular weight scaling method. 

Once the value of Tg has been determined, the values of and i?g for all the molecular weights 
considered here can be evaluated from values of Ru and Rg measured at the various temperatures. 
In the case of Ru, since Tg lies within the range of experimental temperatures, i?^ can be obtained 
by linear interpolation. On the other hand, since Tg lies outside the range of experimental tempera- 
tures for the Ro measurements, i?g can be obtained by linearly extrapolating values of Rg measured 
at a fixed molecular weight and at different temperatures T > Tg, to the value of Tg. Figure |7] 
demonstrates this procedure for two different molecular weights, using values of R^ measured at 
various temperatures. The values of R^ and 7?g determined by this procedure are plotted against 
M in Figure [H which confirms that indeed ideal chain statistics are obeyed at the estimated 6- 
temperature. Estimated values of i?^ and i?g for all the molecular weights considered here, are 
listed in Tables IV and V of the supplementary material. 

Since the scaling of both i?g and i?^ with M at the 0-temperature is the same (M"^), their ratio 
should be a constant. As is well known, experimental observations and theoretical predictions 
indicate that U^^ = RyR^^ is inf act a chemistry inde pendent universal constant (for a recent 



compilation of v alues see Table I in 



1.47934 [IOttingeilll996 



Kroger et al. 



ZimmL 



pOOOll '). Zimm theory predicts a universal value 



195611 . Since we have estimated both R^ and R^, we can 



calculate for all the molecular weights used in this work. The molecular weight independence 
of is shown in Figure |9l It is clear that the values obtained here are close to the value predicted 
by Zimm. 



B. The universal swelling of DNA 

In order to establish the universal nature of the swelling of DNA molecules in solution, it is 
sufficient to show that the data in the crossover region between 6 and athermal solvents can be 
collapsed onto master plots when represented in terms of the scaling variable r. However, it is 
more useful for comparisons between theories and experimental observations for DNA solution 
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FIG. 7. Estimation of /?g by extrapolation of Rg values measured at various temperatures to the value of 
Tq = 13°C, for two different molecular weights. 
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FIG. 8. The variation of the radius of gyration (/?g) and the hydrodynamic radius (/?^) with molecular 
weight (in bp) at the 0-temperature, which is estimated to be 13°C. 
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FIG. 9. The molecular weight independence of U^j^ - The va 



prediction of ?7^„ ~ 1 .47934 in the long-chain limit llQttinger 



ues are in accordance with the Zimm model 



19961 : 



Zimm . 



195611 . 



behaviour, if universal behaviour is established in a framework that enables such a comparison. 
As discussed earlier in section H carrying out a systematic comparison of experimentally deter- 
mined expansion factors with the predictions of two-parameter theory requires relating the scaling 
variable r to the solvent quality parameter z, since theoretical predictions are made in terms of z. 

Many versions of the two-parameter theory can be used for the purpose of comparison with ex- 
perimental data. Early versions of the two-parameter theory, including Flory type fifth power 
the ories, and theorie s based on perturbation expansions, are reviewed in the classic textbook 
by lYamakawal [1200 ill . The use of two- parameter theory to describe numerical results on self- 



avoiding walks on lattices is reviewed in 



Barrett et al. 



[1199 ill . While quantitative agreement be- 



tween experiments and theory for the swelling of static properties suc h as the r adius of gyra 



tion Rg can be obtained with the help of renormalization group theories [|Schafer , 



1 99911 ■ similar 



agreement has not been achieved for the swelling of the dynamic property i?H within this frame- 



work. Recently, quantitative pred ictions of both an c 
Brownian dynamics simulations [iKumar and Prakash , 



aH have been achieved through th e use of 



2003 : 



Sunthar and Prakashl . 



2006ll . These 



simulations, which make use of a narrow Gaussian potential to represent excluded volume inter- 
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FIG. 10. Swelling of Rg vs z- The filled coloured symbols represent experimental data for DNA. BDS 



200311 . with the curve 



refers to the predictions of Brownian dynamics simulations MKumar and Prakashl . 
representing the function /g(z), with constants a = 9.528, b = 19.48, c = 14.92, and m = 0.1339. Black 
squares represent experimental data of po lystyrene in cyclohexane and benzene, reproduced from Miyaki 



and Fujita, 1981 iMiyaki and Fujita , 



198 Ih . 



actions (that enables controlled simulations in terms of the parameter z), give rise to parameter 
free predictions, independent of finite size effects, since they are based on the extrapolation of 
finite chain results to the long chain limit. In particular, the quantitative prediction of was 
show n to require the inclusion of fluctuating hydrodynamic interactions [[Sunthar and Prakashl . 
20061] . Interestingly, it turns out that lattice simulations, predictions of renormalization group 
theory, and results of Brownian dynamics simulations for the swelling of any typical generic prop- 
erty, say a, can be described by universal functions that have the functional form a = f{z), where, 
f{z) = (l+az + b z^ + c ^3-)m/2^ ^-^^^ values of t he constants a, b, c, m, etc. , dependent on the par- 



ticular theory adopted [[Domb and BarrettL 



1976 : 



Kumar and Prakashl . 



2003 : 



Schafed . 



199911 . Here, 



we choose to compare experimental results for ag and an with the results of Brownian dynamics 
simulations, since they represent exact numerical simulations, and are not based on approxima- 
tions unlike earlier analytical theories. The values of the various constants in the functions fg and 
/h that fit the results of Brownian dynamics simulations, are reported in the captions to Figs. \T0\ 
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FIG. 11. Swelling of /?h vs z- The filled coloured symbols 



represent experimental data for DNA. BDS 



refers to the predictions of Brownian dynamics simulations jSunthar and Prakash . 



200611 . with the curve 



representing the function /h(z), with constants a = 9.528, b = 19.48, c - 14.92, and m = 0.0 67. Black stars 



represent several experimental data on synthetic polymers collated in 



Tominaga et al. 



1 200211 ■ 



andm 

The chemistry dependent constant k that relates the value of z in Brownian dynamics simula- 
tions to the experimental value of t remains to be determined. This is done as follows. Consider 
a^g^^ to be the experimental value of swelling at a particular value of temperature T and molecular 
weight M. It is then possible to find the Brownian dynamics value of z that would give rise to 
the same value of swelling, from the expression z = f~^(a^g^^), where /g~' is the inverse of the 
function /g. Since z = kf Vm, where f = (l - ^j, it follows that a plot of /g~^(Qrg''^')/ ^^M versus 
t, obtained by using a number of values of org"'" at various values of T and M, would be a straight 
line with slope k. Once the constant k is determined, both experimental measurements of swelling 
and results of Brownian dynamics simulations can be represented on the same plot. As discussed 
in greater detail in the supplementary material, by using ^-temperature of 13°C for the solvent 
used in this study, the value of k, found by this procedure is 7.4 x 10"^ + 3.0 x 10"'*. It follows that 
for any given molecular weight and temperature, the solvent quality z for the DNA solution can be 
determined. 
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FIG. 12. Comparison of measured values of au against ag. Filled symbols represent experimental data for 
DNA at good solvent conditions (7 > Tg). BDS refers to the predictions of Brownian dynamics simula- 



tions iSunthar and Prakash . 



collated in lTominaga et al. 



120021]. 



2006]. Open circles represent several experimental data on synthetic polymers 



The solvent quality crossover of and qth for DNA obtained by the procedure described above 
are shown in Figs.fTOlandfTTl along wit h Brownian dyna r nics si mulations predictions. Experimen- 



tal data of Miyaki and Fujital [|l98l|] and lTominaga et al\ [|2002|l . which are considered to be highly 



accurate measurements of synthetic polymer swelling, are also plotted in the same figures. An al- 
ternative, z independent, means of presenting the swelling data is displayed in Fig. [121 where the 
values of ag, at particular values of T and M, are plotted versus qth at the same values of T and M. 
From these figures it is evident that, just as in the case of synthetic polymer solutions, irrespective 
of solvent chemistry, the swelling of DNA is universal in the crossover region between theta to 
good solvents. With the experimental swelling of DNA mapped onto theoretical predictions, it 
follows that other instances of universal behaviour, obtained for instance by simulations or theory 
for a particular value of z, can be compared with experimental results for DNA in the solvent used 
here, at the same value of z. 
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FIG. 13. Intrinsic viscosity of 25 kbp and A DNA (48.5 kbp) at 20°C and 25°C, respectively, obtained by 
extrapolating data for In (77o/?/.s)/c and rjpo/crjs at several finite concentrations to the limit of zero concen- 
tration. 

IV. SOLVENT QUALITY CROSSOVER OF THE ZERO SHEAR RATE VISCOSITY 

The intrinsic viscosity [rj] of the various DNA solutions used in this work has been obtained 
from the following two expressions, 

ln(?7o/?7s) 



[rj] = lim ■ 



>o c 



W=lim— (6) 



where, 770 is the zero shear rate solution viscosity, 77, the solvent viscosity, and 77^0 = ijo - Us, the 
polymer contribution to the zero shear rate solution viscosity. Essentially, data for In (77o/?7.v)/c and 
Hpo/cJls obtained at several finite concentrations is extrapolated to the limit of zero concentration, 
as exemplified in Figure [I3] for 25 kbp and A DNA (48.5 kbp) at 20°C and 25°C, respectively. 
From the expression, 77 = 77^ [1 + [77] c + •••], we expect that in the dilute regime (c < c*), 

^ = ^ (7) 
?7,[?7]c* c* 

Figure [13] shows that indeed for all the temperatures measured here, both the 25 kbp and A DNA 
samples exhibit a linear dependence on the scaled concentration c/c* in the dilute regime. 
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FIG. 14. Linear dependence of the viscosity ratio r]po/{r]s[T)]c*) on the scaled concentration c/c*, in the 
dilute solution regime, for 25 kbp and A DNA (48.5 kbp) at a variety of different temperatures. 

According to blob theory, the polymer contribution to the zero shear rate viscosity obeys the 
following scal ing law in the semidilute une ntangled regime in the limits of ^-solvents and very 



good solvents jRubinstein and Colby . 



2003|], 



C \l/(3v-l) 
C* 



(8) 



where, the scaling exponent v is equal to 0.5 in ^-solvents, and 0.588 in very good solvents. The 
semidilute unentangled regime is expected to sp an the range of concentrations from c/c* = 1 to 10 



[|GraessleyLll980l : lRubinstein and ColbyLl2003|l . Figure [T?] displays the ratio ?7po/?7*o (where 77* ^ is 



the value of rjpQ at c = c*) at the ^-temperature estimated here. Using this ratio to represent the 
zero shear rate viscosity data ensures that J]po/ri*po = 1 when c/c* = 1 for all the systems studied 
here. The slope of 2 that is observed in the semidilute regime is consistent with the blob scaling 
law for V = 0.5. T4 DNA, which is the longest molecule in the series studied here, appears to 
follow the blob scaling for the largest concentration range. On the other hand, the 25 kbp and A 
DNA seem to crossover into the entangled regime beyond a concentration c/c* « 3. 

The key result of the present work is displayed in Figure [161 The concentration dependence of 
three different molecular weights of DNA is presented for three different values of the solvent qual- 
ity z. hi order to maintain the same value of solvent quality across the various molecular weights, 
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FIG. 15. Dependence of the viscosity ratio T]po/r]*^Q (where Ty*^ is the value of ripo at c = c*) on the scaled 
concentration c/c*, for 25 kbp, A and T4 DNA, at the ^-temperature. 

it is necessary to carry out the experiments at the appropriate temperature for each molecular 
weight. This procedure would not be possible without the systematic characterisation of solvent 
quality described earlier in section Hill Remarkably, the data appear to collapse onto universal 
curves, independent of DNA molecular weight. Indeed the data appear to satisfy a power law 
analogous to the blob scaling power law. 



77, \c*/ 



l/(3yeff(z)-l) 

(9) 



where, Veff is an effective exponent that is a function of the solvent quality z. The value of Veff (z) can 
be extracted by equating the slope of the (ripo/rj*^^) versus (c/c*) line in Figure [H] to 1/(3 Veff - 1), 
at each value of z. The effective exponent obtained in this manner can be compared to an al- 
ternative effective expon ent that is commonly used to characterise the crossover from 9 to good 



solvents llSchafei] . 



199911 



Although the scaling of Rg with chain length N obeys a power law only in the limit of long 
chains, one can make the ansatz that, at all values of chain length, ~ N"^"^, with Veff representing 
an effective exponent that approaches its critical value as N ^ 00. It is straight forward then to 
show that, 

1 idlnal 

Veif(z) = T + T^r^- (10) 

2 4 omz 
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FIG. 16. Dependence of the viscosity ratio rjpolrf^ on the scaled concentration c/c* in the semidilute 
regime, for 25 kbp, A and T4 DNA, at fixed values of the solvent quality z. 

The curve fitting function for cfg can consequently be suitably differentiated to obtain the z de- 
pendence of the effective exponent Vgff. In an entirely analogous manner, one can make an ansatz 
regarding the scaling behavior of both and with A^. Here, the swelling of the intrinsic vis- 
cosity is defined by the expression, a^, = ([?7]/[?7]e)^''^, where \v1\0 is the intrinsic viscosity at the 
^-temperature. The dependence of Veff on z in each of these cases, for the exact results of Brownian 
dynamics simulations, is illustrated in Figure [iTl The black symbols represent the experimental 
values of Vgff obtained from Figure [16] as described above. It is clear that the experimental Vgff is 
bracketed by the three simulation estimates, with perhaps the closest proximity to the Vgff versus z 
curve obtained from or^. 



Liu et al. 



[1200911 have recently examined the concentration dependence of the longest relaxation 



time Al of semidilute solutions by studying the relaxation of stretched single T4 DNA molecules 
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FIG. 17. Effective exponent Veff obtained from the universal crossover behaviour of ag, a^, and predicted 
by dilute solution Brownian dynamics simulations, compared with the effective exponent extracted from the 
crossover behaviour of the viscosity ratio in the semidilute regim e (displayed in Figure [T6l ). The simulation 



results for are reproduced from a currently unpubhshed work BAhirwal et al. 



201211 



after extension. They find that the scaling of A\ in the semidilute unentangled regime obeys the 
power law expression predicted by blob theory, 



(2-3y)/(3y-l) 



(11) 

It is common to define an alternative larg e scale relaxati on time /l^ based on the zero shear rate 
viscosity rjpQ by the following expression I Ottingei . ll996 1. 



cNAksT 



(12) 



where, A^^ is Avagadro's number and kg is Boltzmanns constant. It is straight forward to show 
that in the semidilute unentangled regime Aj^ obeys the same scaling with concentration as /li. 
As a result, one expects the ratio Aj^l Ai to be independent of concentration. Indeed, in the dilute 
limit it is known that thi s ratio has a universal value, predictions of which range from 1.645 to 
2.39 [ Kroger g/l 2000 1. depending on the theoretical approximation and solvent condition. 
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FIG. 18. Universal ratio of X\ (measured by 



Liu et al. 



1 200911 ) to A„ (measured in the present work) for T4 



DNA at a value of solvent quality z = 2.3, for a range of concentrations spanning the dilute and semidilute 
regime. 



Since Liu et al.\ [|2009l] used the same solvent as the one used in the present study, it is possible 
to estimate the solvent quality of the T4 DNA solution used in their work, which we find to be 



z = 2.3. As a result, by using the value of rjj , 



Liu et al. 



mea sured at the same solvent quality, it is possible 



[I2OO9I1 to Ari measured in the present work at z = 2.3. 



to find the ratio of Ai measured by 
Results for this ratio at a range of concentrations are shown in Figure \TS\ Clearly the ratio has an 
universal value Ar^/Ai «2.369atz = 2.3, independent of concentration. 



Recently. ISomani et all [120 1011 have predicted the dependence of the ratio Aj^/Ai on the solvent 
quality z in the dilute limit with the help of Brownian dynamics simulations. Their results are 
reproduced as the solid curve in Figure \T9\ Since Arj/Ai does not depend on concentration, the 



same curve should be valid at all concentrations. 



estimated here with the results of [|Somani et ali 



t is po ssible therefore to compare the ratio Ar,/Ai 



2OIOI1 . The red symbol in Figure [19] is the value 



of the ratio estimated here at z = 2.3. As can be seen, the estimated value is reasonably close to 
the simulation prediction at this value of z. 
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FIG. 1 9. Dependence of the ratio A,^/Ai on the solvent quahty z in the dilute limit predicted by lSomani et al. 



1120 1011 with the help of Brownian dynamics simulations (solid c urve). The red sy mbol is the value of the 



ratio estimated here at z = 2.3 using the value of Ai measured by 
present work. 



Liu et al. 



1 20091] and A,j measured in the 



V. CONCLUSIONS 



Static and dynamic properties of DNA in a commonly used buffer solution have been charac- 
terized in terms of parameters used in dilute polymer solution theory. Using known scaling rela- 
tionships from polymer solution theory, the 0-temperature (Tg) in this solvent and the crossover of 
the swelling of the two radii in good solvent conditions {T > Tg) were determined. DNA in this 
solvent (10 mM Tris, 1 mM EDTA, 0.5 M NaCl) has = 13 + 2 °C. As with the regular polymers, 
the dilute solution crossover data of DNA can also be collapsed into a unique function of a scaling 
variable z = k(l - Tg/T) -\Im. This confirms the two parameter theory hypothesis for DNA, that 
the static and dynamic properties of DNA in good solvents can be uniquely determined from just 
two parameters: Radius of gyration at Tg (Rp and the solvent quality z. The scaling functions 
also agree within experimental error with the functions for synthetic polymers and model chains 
in good solvents obtained using molecular simulations. 

The concentration dependence of the zero shear rate viscosity of semidilute solutions of un- 
entangled DNA molecules, across a range of molecular weights, is shown to obey the power 

25 



law scaling predicted by blob theory at the ^-temperature. In the crossover region between the 
^-temperature and very good solvents, it is seen that the concentration dependence also obeys 
universal scaling behaviour which can be described by the blob scaling law, but with an effective 
exponent Vgff that depends on the solvent quality in place of the Flory exponent v. 

The results obtained here demonstrate that it is essential to represent observations in terms of 
the solvent quality variable z in order to properly understand and characterise the concentration 
and temperature dependent dynamics of semidilute polymer solutions. 
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Appendix A: Estimation of the persistence length 

At the ^-temperature, since a polymer chain obeys Gaussian statistics, the persistence length P 
is related to the radius of gyration i?g and the contour length L, through the relation P = 3 {Rp^/L. 
In the asymptotic limit (Rg/L) «; 1, this relationship is also valid for a worm- like chain. The 
persistence lengths evaluated using this expression, at the 0-temperature, from known values of 
and L, are shown in Figure [201 It is evident from the figure that the persistence length of DNA at a 
salt concentration of 0.5 M (Na^) is constant across all the molecular weights, with a mean value 
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FIG. 20. Molecular weight independence of DNA persistence length. 
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FIG. 21. Dependence of DNA persistence length on salt concentration, collated from data reported previ- 
ously. The DNA molecular weights used in all these studies range from 3-300 kbp. The mean persistence 
length for all the molecular weights shown in Figure |20] is denoted here by the opaque red circle. Abbre- 
viations: 'FD' - Flow Dichroism, 'FB' - Flow Birefringence, 'MB' - Magnetic Birefringence, 'FE' - Force 
Extension using optical tweezers. 
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that is approximately 44 nm. 

The effects of salt concentration on the persistence length of DNA has been stu died earlier 



through a variety of techniques such as Fl ow Dichroism [Rizzo and Schellmanl 



fringence [|Caimey and Harringt 



X)n[ ll982ll . 



198in. 



Magneti c Birefringence [ IMaret and WeillL 



1997 



'low Bire- 



Wenner et al. 



mm an d 



200211 ■ 



force-extension experiments using optical tweezers llBaumann et ali 
Data from a number of these studies has been collated in Figure [2T1 with the mean value measured 
in the present study included. Clearly, the persistence length in all the earlier studies appears to 
reach an approximately constant value of 45-50 nm for salt concentrations > 0.1 M, suggesting 
that the charges have been fully screened in this concentration regime. Both the value of the per- 
sistence length in the high salt limit, and the threshold concentration for charge screening obtained 
in this work are consistent with earlier observations. 



Some e a rly studies of the ionic depende nce of DNA persistence length 



Borochov et al. 



Kam et al. 



1119811] : 



il98in : lSobel and Harpsti [1199111 . investigated via light scattering, have shown a strong 
dependence of persistence length on salt concentration even up to 4 M. Notably, these measure- 
ments were not at the 0-temperature, and consequently, required both an estimation of the solvent 
quality, and a correction to account for excluded volume effects. The excluded v o lume c orrections 



applie d in the original papers wer e modified in subsequent studies iHagermanl 11198811 : [Manning 



mm-. 



Schurr and AUisonl il98l|] . These aspects might explain the differences observed from 



the present study, where both diffusivity and persistence length measurements suggest that charge 
screening is complete above 10 mM NaCl. 



Appendix B: Thermal blobs and measurements in poor solvents 

The focus of the experimental measurements in the dilute limit reported in the present paper 
is twofold: (i) determining the ^-temperature, and (ii) describing the 6 to good solvent crossover 
behaviour of a solution of double- stranded DNA. The analysis of molecular weight scaling under 
poor solvent conditions has been carried out essentially only in order to locate the 0-temperature. 
As is well known, the experimental observation of single chains in poor solvents is extremely 
difficult because of the problem of aggregation due to interchain attraction. Nevertheless, in this 
section we show that a careful analysis of the data displayed in Figure |3l in the light of the blob 
picture, enables us to discuss the reliability of the measurements that have been carried out here 
under poor solvent conditions. 
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According to the blob picture of dilute polymer solutions, a polymer chain in a good or poor 
solvent can be considered to be a sequence of thermal blobs, where the th ermal blob denotes the 



length scale at which excluded volume interactions become of order k^T [[Rubinstein and Colby . 



200311 ■ Under good solvent conditions, the blobs obeying self-avoiding-walk statistics, while they 



are space filling in poor solven ts. As a result, the mean size R of a polymer chain is given 



by IIRubinstein and Colby . 



200311 ■ 



R = i?biob(r) I I (Bl) 



where, Nk is the number of Kuhn-steps in a chain, //biob is the number of Kuhn-steps in a thermal 
blob, and i?biob is the mean size of a thermal blob. The Flory exponent v is « 0.59 in a good 
solvent, and 1 /3 in a poor solvent. The size of the thermal blob is a function of temperature. For 
instance, under athermal solvent conditions, the entire chain obeys self avoiding walk statistics, 
so the blob size is equal to the size of a single Kuhn-step. On the other hand, for temperatures 
approaching the 0-temperature, the blob size grows to engulf the entire chain. From the definitions 
of Hh and Hg (see Eqs. dU), it is clear that on length scales smaller than the blob length scale 
both these quantities must remain constant. On the other hand, Eq. (|B1|) suggests that on length 
scales large compared to the blob length scale, Hh and Hg must scale as °^ in good solvents, 
and M~^^^ in poor solvents. Figure |22] re-examines the logllH versus logM data for the solvent 
used in this study in the light of these arguments. It is clear that after an initial regime of constant 
values, there is a crossover to the expected scaling laws in both the good and poor solvent regimes. 
The crossover from one scaling regime to the next begins approximately at the blob length scale, 
an estimate of which can be made as follows. 

The requirement that the energy of excluded volume interac tions within a thermal blob a re of 



order ksT leads to the following expressions for A^biob and i?biob [[Rubinstein and Colby . 



2003|], 



^biob(r) =-|^ (B2) 

v{T) 

RuobiT) =— ^ (B3) 
\v(T)\ 



where, bt is the length of a Kuhn-step, and v(T) is the excluded volume at temperature T. The 
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FIG. 22. Different scaling regimes for the scaled variable ITh as a function of molecular weight M. The 
filled red circles correspond to the molecular weight Mbiob of the chain segment within a thermal blob. 

excluded volume can be shown to be related to the temperature through the relation, 

uo |l - ^ j for good solvents, 



v{T) 



(B4) 



for poor solvents. 



where, vq is a chemistry dependent constant. These expressions are consistent with the expec- 
tatio n that v ^ vg in an atheim al solvent {T oo), and v -Vq in a non-solvent (T 



0) [[Rubinstein and Colby . 



200311 . Since measurements of the mean size (via R„ and Ru) have 



been carried out here at various temperatures, and we have estimated both Tg and bt, it is possible 
to calculate v(T) in both good and poor solvents using Eqs. (IBll) to (IB3I) . As a result the size of a 
thermal blob as a function of temperature can also be estimated. 

The equations that govern the dimensionless excluded volume parameter Vo/bl and the molec- 
ular weight Mbiob of a chain segment within a thermal blob, in good and poor solvents, are 
tabulated in Table H when the hydrodynamic radius Ru is used as a measure of chain size. 
Here, m^- is the molar mass of a Kuhn-step, and the universal amplitude ratio Ur has been 



30 



solvent quality 


good 


poor 


^0 

(for M > Mbiob) 


1 

'aUu{URURD)ml~\2v-i \ 


'aUu(URURn)T^ 1 






Mbiob(r) 











TABLE I. Equations for the dimensionless excluded volume parameter vo/b^^, and the molecular weight of 
the chain segment within a thermal blob M^iob, in good and poor solvents. Here, m^- is the molar mass of a 
Kuhn-step, and Ur and Urd are universal amplitude ratios, such that R = URRg, and Kg = UrdRu- 



used to relate the mean size R to Rg (R = URRg), while the universal ratio Urd relates Rg to 
(^g = UrdRh)- The values of these ratios are k nown analytically f or the case of Gaussian 
chains and Zimm hydrodynar nics under 0-conditions IIDoi and EdwardsL 1198611 . and numerically 
in the case of good solvents IIKumar and Prakashl . 1200311 and fluctuating hydrodynamic interac- 



tions [ISunthar and Prakash 



200611 ■ 



Using the known values of a, IIh, bk, nij^, Urd, Ur in the appropriate equations in TableHl we find 
that for sufficiently high molecular weights, vo/bl « 4.7 + 0.5 in both good and poor solvents. This 
is significant since an inaccurate measurement of mean size in a poor solvent (as a consequence 
of, for instance, chain aggregation), would result in different values of Vo/bl in good and poor 
solvents. Further evidence regarding the reliability of poor solvent measurements can be obtained 
by calculating M^,lo^,{T) in good and poor solvents. 

Figure [23] displays the variation of Mbiob with respect to the temperature difference T - Ty, cal- 
culated using the equations given in Table IB The figure graphically demonstrates the temperature 
dependence of the blob size, and confirms that essentially the blob size is the same in either a 
good or poor solvent when the temperature is equidistant from the ^-temperature. The symbols in 
Fig. [23] denote values of Mbiob, evaluated at the temperatures at which experimental measurements 
have been made. These values have been represented by the filled red circles in Fig. [22] As can be 
seen from Fig. [221 the magnitude of Mbiob is roughly consistent with the location of the crossover 
from the scaling regime within a blob, to the scaling regime that holds at length scales larger 
than the blob, in both good and poor solvents. The two scaling regimes, at 25°C, are illustrated 
explicitly in Fig. [24] 
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FIG. 23. Variation of the molecular weight of the chain segment within a thermal blob with respect to tem- 
perature, on either side of the ^-temperature (see Table U for the equations governing Mbiob)- The symbols 
denote values at temperatures at which experimental measurements have been made. 
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FIG. 24. The variation of hydrodynamic radius (/?h) with molecular weight (in bp) at 25°C. 
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FIG. 25. Spinodal curves and critical temperatures and concentrations (filled circles) predicted by Flory- 
Huggins mean-field theory for a range of molecular weights. Values of concentrations and temperatures at 
which the poor solvent experiments have been conducted are also indicated. 



The possibility of phase separation under poor solvent conditions, as polymer-solvent inter- 
actions become less favourable, is the primary reason for the difficulty of accurately measuring 
the size scaling of single chains. An approximate estimate of the thermodynamic driving force 
for phase separation can be obtained with the help of Flory-Huggins mean field theory. Since 



the F 



ory-Huggins y pararneter is related to the excluded volume parameter through the rela- 



tion [[Rubinstein and Colby . 



2mx 



^ |l - ^j, and we have estimated the value of wq/^^ 
in both solvents, the phase diagram predicted by Flory-Huggins theory for dilute DNA solutions 
considered here can be obtained. It is appropriate to note that we are not interested in accurately 
mapping out the phase diagram for DNA solutions with the help of Flory-Huggins theory. This 
has already been studied in grea t detail, using sophistic ated versions of mean-field theory, start- 



ing with the pioneering wor k of 



is an active field of research [|Teif and Bohin 



linJ . 



Post and Zimml 1119821]. and the prob 



2011 



Yoshikawa et al. 



em of DNA condensation 



1996 



201 1|1 . Our primary 



interest is to obtain an approximate estimate of the location of the current experimental measure- 
ments relative to the unstable two-phase region (whose boundary is determined by the spinodal 
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curve), since phase separation can occur spontaneously within this region. Figure [25] displays the 
spinodal curves for the 25 to 289 kbp molecular weight samples, predicted by Flory-Huggins the 



ory. D etails of how these curves can be obtained are given, for instance, in [Rubinstein and Colby 



ll2003ll . Also indicated on each curve are the critical concentration and temperature. It is clear 
by considering the location of the symbols denoting the concentration-temperature coordinates of 
the poor solvent experiments, that for each molecular weight, they are located outside the unstable 
two-phase region, lending some justification to the reliability of the present poor solvent measure- 
ments. It is appropriate to note here that mean-field theories do not accurately predict the shape 
of the binodal curve, and in general concentratio n fluctuations tend to make the curve wider close 
to the critical point [Rubinstein and Colbyl . 1200311 . Interestingly, even for the 289 kbp sample, that 
has a very large molecular weight (« 1 .9 x 10^ Dalton), there is still a considerable gap between the 
critical and ^-temperatures (« 4°C). The reason for this is because the stiff"ness of double- stranded 
DNA leads to a relatively small number of Kuhn-steps (1 120) even at this large value of molecular 
weight, and the value of the critical temperature predicted by Flory-Huggins theory depends on 
the number of Kuhn-steps in a chain rather than the molecular weight. 
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1 



I. DETAILS OF STRAINS, WORKING CONDITIONS AND PROCEDURES FOR 
PREPARING LINEAR DNA FRAGMENTS. 



Recently, a range of special DNA c onstructs, from 3 - 300 kbp, have been genetically engi 



neered into bacterial strains of E. coli ULaib et a/.L 1200611 . which can be selectively extracted for 
rheological studies. Primarily they fall into three categories: plasmids, fosmids and Bacterial Ar- 
tificial Chromosomes (BAC). Altogether si x samples (two p lasmids, two fosmids and two BACs), 



which were originally prepared elsewhere IILaib et alx 1200611 were procured from Dr. Brad Olsen, 



the three 



20060. In 



California Institute of Technology, USA. Throughout the work, the nomen clature of a l _ 
types of DNA samples has been used as in the originally published work [ILaib et alx 
addition, two special bacterial strains containing the plasmids: pBSKS (2.7 kbp) and pHCMCOS 
(8.3 kbp) were provided by Dr. S. Noronha, Dept. of Chemical Engineering, IIT Bombay (In- 
dia). The details about size, growth conditions of bacteria and single cutters of the DNA samples 
are mentioned in Table HI After procurement of samples (in the form of agar stab cultures of E. 
coli), glycerol freeze stocks were made using 50% glycerol and stored at -80°C. The cultures 
can be stored in th is way for several years and can be used at any time to produce DNA samples 



[ILaib et al. 



200611 



Standard procedures [[Laib et al. 



2006 



Sambrook and Russell! 



200111 involving alkaline lysis 



(mediated by NaOH) were adopted for extraction, linearization and purification of plasmids, fos- 
mids and BAC from the cultures. For high copy number plasmids, no inducer was added. For low 
(fosmids) and very low (BACs) copy number samples, L-arabinose was added as inducer. From 
each freeze stock, 15/il of ice was scrapped and transferred to 40 ml LB medium with proper an- 
tibiotic (as mentioned in TableH]) and incubated overnight (16 - 18 hours) at 37°C with vigorous 
shaking (200 - 250 rpm). The overnight grown culture was poured into microcentrifuge tubes and 
cells were harvested by centrifugation. The bacterial pellet (obtained above) was resuspended in 
100/zl of ice-cold Sol ution I (4°C) followed by 200^ 1 of freshly prepared Solution II and ISQjA of 



ice-cold Solution III [Sambrook and Russell! 



200 ill . The tubes were stored on ice for 3-5 min- 



utes and centrifuged. The supernatant was transferred to a fresh tube. The precipitate (containing 
mainly the cell debris and genomic DNA) was discarded. RNase was added (at 10 //g/ml) to the 
tube and incubated at 37°C for 20 minutes. Equal volume of Phenol-Chloroform-Isoamyl Alcohol 
(25:24:1) mixture was added and mixed well by vortexing. After centrifugation, supernatant was 
transferred to a fresh tube. Equal volume of Chloroform was added and centrifuged. The super- 



TABLE I. DNA Fragments. Here 'LB' stands for Luria Bertini broth, 'Ant refers to Antibiotic resistance, 
'Amp.' refers to Ampicillin, 'CAM' refers to Chloramphenicol, 'Kan' refers to Kanamycin, all the cul- 
tures were incubated overnight at 37°C with vigorous shaking (200-250 rpm). L-arabinose (inducer) was 
used at a concentration of O.Olg per 100 ml (stock concentration: 5g in 100ml). Stock concentrations for 
preparations of Ampicillin, Chloramphenicol and Kanamycin were 100 mg/ml, 25 mg/ml and 100 mg/ml 
respectively. The working concentrations for Amp., CAM and Kan are 100 //g/ml, 12.5jUg/ml and 100//g/ml 
respectively. Growth conditions for all the plasmids are same (LB + Amp.) except pHCMC05 (LB + Amp. 
+ CAM). For both the fosmids, growth conditions are identical (LB + CAM + L-arabinose). For both the 
BACs, growth conditions are same (LB -i- CAM -i- Kan -i- L-arabinose). 



Type 


Name 


Size (kb) /(Notation) 


Ant'^ 


1 Cutter 


Plasmid 


pBSKS(-F) 


2.9/F2.9 


Amp. 


BamHI 




pYES2 


5.9/F5.9 


Amp. 


BamHI 




pHCMC05 


8.3/F8.3 


Amp. + CAM 


BamHI 




pPIC9K<TRL5> 


11.1/Fll.l 


Amp. 




Fosmid 


pCClFOS-25 


25/F25 


CAM 


Apal 




pCClFOS-45 


45 / F45 


CAM 


Apal 


BAC 


CTD-2342K16 


114.8/F114.8 


CAM + Kan 


Mlul 




CTD-2657L24 


289 / F289 


CAM + Kan 


Mlul 



natant was transferred to a fresh tube. Two volumes of chilled 100% ethanol (at 4°C) was added at 
room temperature kept for 7-8 hours at -20° C. The tube was then centrifuged and the supernatant 
removed by gentle aspiration. The tube was kept in an inverted position in a paper towel to allow 
all of the fluid to drain away. Following this, 1 ml of 70% ethanol was added to the tube and 
centrifuged. When all of the ethanol was evaporated, the resulting DNA pellet was dissolved in 
50yul of Milli-Q grade water and stored at -20°C. 

To linearize the extracted DNA fragments, 39/il of water was added to a 1 .7 ml microcentrifuge 
tube, followed by lO/A of corresponding lOX Assay Buffer (working concentration is IX) and 
50fA of DNA solution (purified DNA stored at 4°C). 1 lA of appropriate enz yme was added. A 



thumb rule is 0.5 - 1 U enzyme for 1 mg DNA [ISambrook and Russellll200lll . The samples were 



mixed well with micropippette (wide bore tips) for several times. The reaction mix (100//1) was 



TABLE II. Composition of the solvent used in the study. 
10 mM Tris 
1 mM EDTA 
0.5 M NaCl 
Water 



incubated at 37°C for three hours. After restriction digestion / linearization, it is necessary to 
remove the enzymes / other reagents present in the reaction mix so that they do not interfere with 
the downstream application/s like light scattering studies, rheometry etc. For this normal phenol- 
chloroform extractimJo]J^wed_by^^ precipitation of DNA was carried out as described 



elsewhere [ Sambrook and Russell . 



2001 1. 



II. SOLVENT COMPOSITION 

The solvent used in this study (see Table HH) is primarily the widely used Tris-EDTA buffer, 
supplemented with 0.5 M NaCl. The measured viscosity of the solvent at 20°C is 1.01 mPa-s. 



III. QUANTIFICATION OF DNA SAMPLES 



After the DNA samples were extracted and purified, their purity were determined using the 
Nano-Photometer (UV-VIS Spectrophotometer, IMPLEN, Germany). Optical Density (O.D.) 
readings were taken at three different wavelengths: 260 nm, 280 nm and 23 nm. The ratio of ab- 
sorbance at 260 nm to that of 280 nm gives a rough indication of DNA purity [[Sambrook and RussellL 
200 in . The concentrations were calculated from abs orbance reading at 260 nm (DN A shows ab- 
sorption peak at 260 nm) by Beer-Lambert's Law llSambrook and Russelll . 1200111 and also by 
agarose gel electr ophoresis through a serial dilution of DNA samples as suggested elsewhere 



jLaib et a/.L 



20061] . All the linear DNA samples demonstrated A260M280 ratio of 1.8 and above 



This indicates good purity for DNA samples, though it is largely an assumption [[Laib et ali 



2006h 



and A2m/-A 2- x() rati o from 2.0 to 2.2 (absence of organic reagents like phenol, chloroform etc) 



[[Laib et al. 



200611 . The low molecular weight linear DNA fragments (plasmids) were quanti- 



fied through agarose gel electrophoresis with a known standard 1 kbp DNA marker (Fermentas). 



TABLE III. c* values (in yUg/ml) for different DNA fragments used in this study at different temperatures. 



Size (kbp) 15°C 20°C 25°C 30°C 35°C 



2.9 


353 


323 


295 


241 


210 


5.9 


244 


231 


210 


154 


141 


8.3 


198 


187 


163 


124 


104 


11.1 


172 


152 


129 


103 


84 


25 


110 


90 


75 


60 


53 


45 


68 


51 


39 


34 


28 


114.8 


31 


22 


15 


12 


10 


289 


16 


8 


5 


4 


3 



For low copy number fragments (fosmids) and very low copy number samples (BACs), it was 
confirmed that the samples were not sheared during extraction by running a very low concentra- 
tion agarose gel for extended period at low voltage. A loss of 25% - 50% was observed in the 
amount of DNA samples after th e linearization procedure. Thi s is attributed to purification steps 



by phenol-chloroform extraction llSambrook and Russell 



200 in . 



IV. SAMPLE PREPARATION FOR LIGHT SCATTERING 



All the linear DNA fragments were dissolved individually in the solvent and were charac- 
terized by dynamic light scattering (DLS) for the hydrodynamic radius, Ru and by static light 
scattering (SLS) for the radius of gyration, R^. For both dynamic and static light scattering, an 



extensive sample preparation method was followed t o ensure repeatability. The methodo 



sample preparatiq i i was modified from earlier studies IILewis et al. 



Sorlie and Pecoral . 



1985 : 



Selis and Pecora . 



ogy of 



1995 



199011 and was repeated before each measurement. The cuvette was washed 



with ethanol (0.5 ml) for 5 times and kept for 15 minutes inside laminar air flow. It was followed 
by wash with milliQ grade water for 10-15 minutes continuously. In the meantime, the solvents 
were filtered with 0.45jU membrane-filter (PALL Corp., USA) with 2 different membranes consec- 
utively. After filtration, DNA was added to make final concentration of c/c* =0.1 (for DLS) and 
c/c* = 0.05 - 0.1 (for SLS). The c* values for different DNA fragments used in this study at 
different temperatures are summarized in Table Hill 



0.7 

A 

+ 0.5 

0.4 

V 0.3 
II 

^ 0.2 

O 0.1 

0.0 

10° lo' lo' lo' id' lO' lO' lO' 
T(HS) 

FIG. 1. Intensity autocorrelation spectra for 25 kbp DNA at 25°C and c/c* =0.1. 

V. CHARACTERIZATION OF LINEAR DNA SAMPLES 

Dynamic Light Scattering was used for measuring the hydrodynamic radii of all the molecular 
weights at different temperatures: 5°C, 10°C, 15°C, 20°C, 25°C, 30°C and 35°C. Static Light 
Scattering was used for measuring the radii of gyration of all the molecular weights at different 
temperatures: 15°C, 20°C, 25°C, 30°C and 35°C. Static Light Scattering was also used to find out 
the second virial coefficients at different temperatures: 11.2°C, 13°C, 14°C, 15°C and 20°C. 
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VI. DYNAMIC LIGHT SCATTERING: TAKING MEASUREMENTS 

The dynamic size characterization of dilute DNA solutions was done using a Zetasizer Nano 
ZS (ZEN3600, MALVERN, U.K.) particle size analyzer with temperature control fitted with a 
633 nm He-Ne laser using back- scattering detection. This instrument uses dynamic light scatter- 
ing (or quasi-elastic light scattering / photon correlation spectroscopy) to measure the diffusion 
coefficient which is then converted to an average hydrodynamic size of particles in solution us- 
ing the Stokes-Einstein equation. A Standard Operating Procedure (SOP) was created using the 
Dispersion Technology Software (DTS 5.00, MALVERN, U.K.) to achieve the desired outcome 
(Ru) without manual intervention. Scattering of the DNA solutions was measured at a fixed 173° 
scattering angle (this enables measurements even at high sample concentrations and the effect of 
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FIG. 2. Application of the Guinier approximation. Intensity as a function of scattering wave vector, mea- 
sured for 25 kbp DNA at 25°C and c/c* = 0.1. 

dust is greatly reduced). The temperature range investigated was from 5 to 35°C. A typical ex- 
ample of measured intensity autocorrelation spectra for 25 kbp DNA at 25°C and c/c* = 0.1 is 
shown in Figure [B Note that, G(t) =< I(t)I(t + r) >, where / is the intensity of scattered light, 
and T is the time difference of the correlator. The Zetasizer Nano ZS has the ability to measure a 
wide size range (0.6-6000 nm in diameter). In this paper, we have reported sizes roughly in the 
range 140-2800 nm in diameter, which is within the size range of the instrument. Readings were 
taken in three temperature scans (a sequence of High-Low, Low-High, and High-Low temperature 
settings); with 5 readings at each temperature. The mean of 15 readings was taken as final data 
point at each temperature for each DNA fragment. 

Measured values of Ru are reported in Table |lVl As expected, the hydrodynamic radius of 
DNA increases with both temperature and molecular weight. The values of i?^ obtained by linear 
interpolation are indicated in italics. 



VII. STATIC LIGHT SCATTERING: TAKING MEASUREMENTS 



The static light scattering (SLS) measurements were obtained from a BI-200SM Goniometer 
(Brookhaven Instruments Corporation, USA) using BI-SLSW Static Light Scattering Software. 
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TABLE rV. Hydrodynamic Radius (Ru) of linear DNA at different temperatures. Each data point corre- 
sponds to the intensity peaks from DLS measurements. The mean of 15 readings was taken as final data 
point at each temperature for each DNA fragment. The values of obtained by linear interpolation are 
indicated in itaUcs. 



Sequence length 


2.9 kbp 


5.9 kbp 


8.3 kbp 


11.1 kbp 


Temperature 


/?H (in nm) R 


H (in nm) /?h (in nm) /?h (in nm) 


5°C 


73+4 


104±3 


123±3 


141±3 


10°C 


77±3 


109±3 


131±3 


152±3 


13°C 


79±2 


115±4 


137±3 


157±3 


15°C 


85±3 


121±3 


145±3 


167±3 


20°C 


87±3 


124±3 


148±3 


173±4 


25°C 


90±3 


131±5 


155±2 


183±6 


30°C 


96±2 


136±4 


162±3 


189±3 


35°C 


101±4 


145±7 


174±3 


203±5 




Sequence length 


25 kbp 


45 kbp 


114.8 kbp 


289 kbp 


Temperature 


/?H (in nm) Ru (in nm) Ru (in nm) /?h (in nm) 


5°C 


203±4 


258±5 


385±13 


540±35 


10°C 


226±5 


303±6 


473±14 


718±46 


13°C 


239±7 


315±9 


503±13 


789±40 


15°C 


258±3 


349±4 


560±18 


897±57 


20°C 


267±8 


367±4 


607±13 


1025±39 


25°C 


286±5 


397±5 


677±15 


1201 ±49 


30°C 


297±4 


417±6 


722±13 


1300±38 


35°C 


313±8 


431 ±8 


753±19 


1363±57 



A separate temperature control system (PolySc, USA) was arranged. For measurements, the 
readings for each molecular weight were taken at different angles at different temperatures. The 
Guinier approximation for calculating was employed. The angle range was selected based on 
sample concentration. A typical example of the application of the Guinier approximation for 25 
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TABLE V. Radius of gyration (7?g) of linear DNA fragments at different temperatures. Each data corre- 
sponds to an average of 10 measurements performed in two temperature scans. The Rg values at Ty = 13°C 
have been estimated by linear extrapolation of data points, and indicated in italics. 



Sequence length 


2.9 kbp 


5.9 kbp 


8.3 kbp 


11.1 kbp 


Temperature 


/?g (in nm) Rg (in nm) Rg (in nm) Rg (in nm) 


13°C 


125±2 


174±4 


206±5 


234±9 


15°C 


129±2 


184±5 


221 ±7 


255±6 


20° C 


133±5 


187±8 


225±8 


266±9 


25°C 


137±6 


193±5 


236+4 


281±8 


30°C 


147±3 


214±4 


258±2 


303±7 


35°C 


154±3 


220±5 


274±4 


324±9 




Sequence length 


25 kbp 


45 kbp 


114.8 kbp 


289 kbp 


Temperature 


Rg (in nm) Rg (in nm) Rg (in nm) Rg (in nm) 


13°C 


342±11 


458±10 


744±13 


1212±37 


15°C 


389±6 


554±15 


980±22 


1667±46 


20° C 


415±5 


612±14 


1098±24 


2102±30 


25°C 


441±5 


666±10 


1242±23 


2426±38 


30° C 


476±7 


696±10 


1347±17 


2664±43 


35°C 


495 ±7 


743±21 


1415±32 


2782±47 



kbp DNA at 25°C and c/c* = 0.1 is shown in Figure [21 The temperature range investigated was 
from 15 -35°C. Readings were taken in two temperature scans; with 5 repeats at each temperature. 
The mean of 10 repeats was taken as final data point at each temperature for each DNA fragment. 
For determining the transition in the second virial coefficient (at the theta temperature), the inten- 
sity of scattered light was measured as a function of scattering vector (which is a function of the 
scattering angle) and DNA concentrations. Rg of the DNA increase with both temperature and 
molecular weight. This is reported in Table IVl The values of i?g obtained by linear extrapolation 
are indicated in italics. 
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-8.4 

-8.7 



c/c*= 0.35 
c/c* = l).4 



9 X 10 



FIG. 3. Variation of reciprocal intensity for 25 kb linear DNA as a function of (at 20°C). Here, IS a 
function of the scattering angle. The readings are extrapolated to find the intensity at zero scattering angle, 

/o. 

VIII. DETERMINATION OF SECOND VIRIAL COEFFCIENTS BY SLS 



In this work, a single, linear, medium molecular weight DNA fragment (25 kbp) is taken and 
the intensity of scattered light, is determined as a function of scattering vector, and poly- 
mer concentration, c through SLS, at temperatures above and below the estimated 0-point. The 
readings from a typical plot of reciprocal scattering intensity (l/Ie) against q^ (see Figure [3]), are 
extrapolated to find the intensity at zero scattering angle (/„). Intensity readings are taken at 5 
different temperatures between 10 - 20°C. Four different concentrations are used: c/c* = 0.2, 0.3, 
0.35 and 0.4 (c* is the overlap concentration, based on our measur ed R^). The second virial coef- 



ficients are determined according to a procedure described in Ref. Qrofino and Mickey Jr. 



■ 11963 



The dependence of red uced scattering ratios at zero angle (c//o) on the reduced concentration 



(c/c*) can be written as jOrofino and Mickey Jr.L 



19631] 



(c/hf^ = (1/M)[1 +A2(c/c*) + A3(c/c*)' + ...] (1) 

Here, M is Molecular Weight, A2 is second virial coefficient and A3 and subsequent coefficients 
resemble interaction terms of higher order. From the plots of (c//o)'^^ against c/c* (see FigureS]), 
slopes (A2/M) are calculated at each temperature in accordance with Equation ([T]) by fitting the 
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FIG. 4. Variation of reciprocal reduced intensity vs concentration . Slopes are calculated at each temperature 
by linear fitting of the data points. 
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FIG. 5. Determination of the chemistry dependent constant k. The data p oints (except 1 14.8 and 289 kbp) 



are least square fitted with a straight line and the slope of this line gives k BKumar and Prakash , 



200311 ■ 



data points with straight lines and neglecting A3 and the other higher order coefficients. 



IX. DETERMINATION OF THE CHEMISTRY DEPENDENT CONSTANT k 



Th e k value for the solvent has been determined by adopting a procedure elaborated in an earlier 
work [IKumar and Prakashl . 1200311 . and briefly summarized in the main text (see Section III B). 

11 



Figure [5] displays plots of /g'^Q^g"'")/ Vm versus f for the solvent used in this study. Only the 
temperatures above the theta point are considered here. The data points were least square fitted 
with a straight line, and the slope k determined, ^-temperature value of 13°C has been used. 
The value of k found by this procedure is 7.4 x lO""* + 3.0 x 10""^. However, while estimating 
the value of k, very high molecular weight linear DNA fragments (114.8 and 289 kbp) were not 
considered, since they deviated from the general behaviour. The reason for this could be that at 
high temperatures (30°C and 35°C), the high molecular weight linear DNA chains swell profusely 
and may start overlapping with each other (in a crossover region between dilute and semi-dilute 
solutions) resulting in solutions that may not be dilute in a physical sense. 
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